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The  pumpj  et  is  a  unique  fluid  machine  which/Mrtil'ffcefi  retarded  wake  flow  and 
produces  high  propulsive  efficiency  such  as  90%.  The  existing  pumpjet  design 
method  is  based  on  a  simple  two-dimensional  graphic  method  which  was  used  for 
pump  design.  As  the  demand  for  the  speed  of  underwater  vehicles, increased  in 
recent  years,  the  existing  design  method  became  inappropriate.  Effort  has  been 
made  to  develop  a  new  three-dimensional  pump  design  method  by  combining  a  blade- 
through  flow  theory  with  blade- to-blade  flow  theory.  Such  a  method  requires 
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20.  many  supporting  sub-theories  to  be  developed.  The  foundation  work  for  the 
blade-through  flow  theory,  i.e.,  streamline  curvature  method,  three- 
dimensional  flow  mapping  technique,  as  well  as  two-dimensional  cascade 
theory,  ^»as  'beefl>  established  in  the  FY-85  GHR/program.  During  this  FY-86 
GHR  program,  the  blade-to-blade  flow  theory  with  corrections  due  to  the 
three-dimensionality  has  been  established.  Particularly,  the  flow  skewness 
which  escaped  many  researchers  has  been  found  essential  for  accurate  design 
of  pumpjet.  ,  ‘ 
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1.0 


BACKGROUND 


The  pumpjet  is  considered  to  be  one  of  the  most  promising 
candidate  propulsors  for  high  speed  underwater  vehicles  and, 
as  a  matter  of  fact,  it  has  recently  been  employed  for  MK-48 
torpedoes,  ALWT — Advanced  Light  Weight  Torpedo  (now  called 
MK-50 ) ,  and  other  underwater  vehicles.  The  pumpjet 
superiority  over  other  propulsion  devices  is  represented  by 
two  major  factors,  i.e.,  high  efficiency  and  quietness. 

The  pumpjet  is  one  of  few  fluid  devices  which  positively 
utilizes  retarded  wake  flow  and  produces  high  propulsive 
efficiency.  This  peculiar  situation  may  be  understood 
readily  by  considering  the  momentum  equati/i  applied  to  a 
control  volume  surrounding  an  underwater  vehicle,  fixed  to 
the  inertial  coordinate  system.  For  a  conventional  pro¬ 
peller,  the  velocity  of  flow  coming  into  a  propeller  blade 
is  approximately  equal  to  the  vehicle  speed  since  the  pro¬ 
peller  diameter  is  large  enough  to  enjoy  the  free  stream 
flow.  In  order  for  the  propeller  to  generate  an  effective 
thrust,  it  should  accelerate  the  flow,  the  ejected  flow 
speed  being  faster  than  the  incoming  flow.  It  means  that  a 
certain  amount  of  the  energy  imparted  on  the  fluid  by  the 
thruster  is  dumped  in  the  surrounding  water.  On  the  other 
hand,  for  a  pumpjet,  the  incoming  flow  velocity  is  retarded 
and  slower  than  the  free  stream  velocity.  In  order  to 
generate  a  thrust,  again  this  flow  has  to  be  accelerated. 
However,  if  the  pumpjet  is  properly  designed,  the  acceler¬ 
ated  flow  velocity  can  be  almost  the  same  as  the  vehicle 
speed.  The  ejected  flow  out  of  the  pumpjet  has  little  rela¬ 
tive  velocity  and  thus  leaves  hardly  any  jet  wake  behind  the 
vehicle.  Compared  to  a  vessel  with  a  conventional  pro¬ 
peller,  a  vehicle  with  a  pumpjet  generates  much  less  wasted 
energy  in  the  flow  field.  This  is  the  major  reason  why  the 
pumpjet  can  produce  a  high  propulsive  efficiency  such  as  90* 
or  higher  when  it  is  properly  designed. 
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Quietness  is  a  guaranteed  aspect  with  the  pumpjet.  As  can 
be  seen  from  its  configuration  (Figure  1-1},  a  long  shroud 
completely  surrounding  the  rotor  helps  prevent  rotor  noise 
from  emitting  into  the  outside  flow  field.  Furthermore, 
this  "internal"  flow  machine  has  better  resistance  charac¬ 
teristics  against  cavitation,  resulting  in  quieter  shallow 
water  operation  where  propulsors  are  most  susceptible  to 
cavitation . 

However,  in  order  to  achieve  such  a  high  standard  of  perfor¬ 
mance,  there  are  many  penalties  to  be  paid  in  reality.  The 
first  such  penalty  naturally  stems  from  the  pumpjet’s  util¬ 
izing  the  velocity-retarded  wake  flow.  A  typical  meridional 
flow  distribution  at  the  inlet  of  pumpjet  rotor  is  shown  in 
Figure  1-2;  the  velocity  at  the  hub  is  only  30%  of  the  free 
stream  velocity  and  rapidly  increases  to  75%  at  the  shroud 
internal  boundary.  This  large  velocity  gradient  in  the 
transverse  direction  is,  of  course,  built  up  by  the  viscous 
boundary  layer  effect  and  is  one  of  the  key  features  causing 
difficulties  in  design,  fabrication  and  eventually  in 
achieving  the  pumpjet  high  performance. 

When  one  designs  an  axial  or  a  near  axial  pump,  it  is  cus¬ 
tomary  to  distribute  the  blade  loading  from  hub  to  tip  in  a 
forced  vortex  or  a  free  vortex  distribution  method,  such  as 
shown  in  Figure  1-3.  Such  distribution  methods  are  impor¬ 
tant  in  obtaining  a  nearly  uniform  discharge  jet  behind  the 
rotor  in  order  to  minimize  the  mixing  loss.  However,  a 
serious  problem  arises  in  attempting  to  implement  either 
forced  vortex  or  free  vortex  loading  distribution  against 
the  flow  field  having  a  large  velocity  gradient  as  shown 
in  Figure  1-2.  Due  to  the  lack  of  enough  meridional  flow 
velocity  near  the  hub,  the  blade  there  should  be  designed 
to  have  extremely  large  incidence  angle  as  well  as  large 
camber.  It  is  for  this  reason  that  the  pumpjet  rotor 
designed  to  date  has  a  distorted  profile  shape  from  hub  to 
tip  (Figure  1-4).  If  this  were  a  conventional  propeller, 


the  stagger  angle  would  become  smaller  towards  the  hub  and 
the  camber  would  stay  more  or  less  constant.  However,  for 
the  reason  mentioned  above,  the  pumpjet  blade  stagger  angle 
first  becomes  smaller  up  to  the  midspan  area  but  becomes 
larger  toward  the  hub  and  thus  the  camber  is  designed  to  be 
substantially  larger. 

This  unusual  rotor  blade  setup  causes  various  hydrodynamic 
problems.  Since  a  typical  flow  incidence  angle  near  the  hub 
has  to  be  extremely  high  (e.g.,  30*),  even  a  slight  error  in 
design  may  cause  flow  separation,  possibly  cavitation  and 
then  noise  generation.  Furthermore,  even  if  design  is  made 
properly,  the  same  vulnerable  situation  is  generated  with  a 
slight  flow  disturbance  or  blade  deformation  due  to  fabrica¬ 
tion  inaccuracy1.  The  existing  one-dimensional  graphic 
pumpjet  design  method  with  empirical  corrections  for  the 
cascade  effect  fails  to  design  a  pumpjet  free  of  flow 
separation.  The  three-dimensional  effect  due  to  the  diag¬ 
onal  flow  configuration  and  the  cascade  effect  are  not  prop¬ 
erly  taken  into  account.  It  is  for  this  reason  that 
development  of  a  more  accurate  three-dimensional  pumpjet 
design  theory  is  in  order. 


1  Some  pumpjet  rotors  are  produced  by  investment  casting 
process  so  that  the  fabrication  accuracy  cannot  be 

expected  to  be  high. 


2.0 


OBJECTIVES 


The  objectives  of  the  work  under  the  GHR  program  are: 

1)  to  develop  a  more  reliable  and  accurate  pumpjet 
design  theory  based  on  a  three-dimensional  pump 
design  concept, 

2)  to  accurately  incorporate  the  cascade  effect  into 
the  theory  and  then, 

3)  to  improve  the  pumpjet  performance  characteristics. 


The  objective  of  the  FY-86  study  was  to  develop  a  mathemati¬ 
cal  model  for  the  blade-to-blade  flow  for  the  three- 
dimensional  pumpjet  design  method  selected  in  FY-85.  In 
FY-86,  the  blade-through  flow  theory  will  be  incorporated 
into  the  blade-to-blade  theory,  resulting  in  forming  a 
complete  three-dimensional  pumpjet  design  theory. 
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THREE-DIMENSIONAL  DESIGN  METHOD 


Design  of  a  pumpjet  for  an  underwater  vehicle  requires  pre¬ 
liminary  information  on  the  vehicle  including  its  geometry 
and  hydrodynamic  drag  coefficient.  .  Furthermore,  most 
importantly,  the  velocity  profile  at  an  upstream  reference 
section  should  be  obtained  either  analytically  or 
experimentally.  Any  error  in  the  velocity  profile  would 
result  in  a  pumpjet  of  lower  efficiency  or  failure  of  the 
pumpjet  meeting  the  specif ications  at  the  design  point.  In 
the  present  study,  it  is  assumed  that  this  velocity  profile 
is  given  at  a  goal  speed  or  at  the  corresponding  Reynolds 
number . 

The  first  step  for  design  of  a  pumpjet  (Figure  3-1)  is  to 
determine  the  shroud  intake  diameter.  From  the  viewpoint 
of  cavitation,  the  maximum  and  minimum  shroud  diameter  to 
prevent  cavitation  must  exist.  If  it  is  too  large,  the 
rotor  blade  tip  speed  becomes  too  high  so  that  cavitation 
occurs.  On  the  other  hand,  if  it  is  too  small,  the  rotation 
speed  must  be  increased  to  generate  the  required  head  so 
that  the  chance  of  cavitation  inception  also  increases. 
Another  aspect  of  determining  the  shroud  diameter  stems  from 
the  consideration  of  overall  propulsive  efficiency.  The 
equation  for  global  momentum  balance  should  be  able  to 
determine  an  efficiency-optimum  shroud  diameter  for  the 
given  velocity  profile  and  vehicle  drag. 

Once  the  shroud  diameter  is  determined,  streamlines  will 
be  calculated  by  using  the  streamline  curvature  method 
C  SC  M  j .  In  this  calculation,  the  loading  distribution  on  the 
rotor  and  blade  thickness  must  be  assumed  in  advance.  One 
of  the  major  concerns  in  using  the  existing  streamline  cur¬ 
vature  method  lies  in  the  fact  that  SCM  may  only  be  used  for 
relatively  uniform  incoming  flow,  but  may  generate  a  sub¬ 
stantial  error  for  a  thick  wake  flow,  i.e.,  highly  retarded 
flow  due  to  the  viscous  boundary  layer  on  the  vehicle  hull. 


Detailed  mathematical  formulation  and  sample  calculations 
have  been  presented  in  the  1985  report  (Furuya  and  Chiang, 
1986).  Also  included  are  discussions  regarding  the  problems 
of  application  of  conventional  SCM  to  the  thick  wake  flow. 


The  next  step  of  the  design  method  is  to  map  the  stream  tube 
or  surface  calculated  by  SCM  onto  a  plane  so  that  the  rotor 
blades  are  mapped  into  cascade  configuration.  If  the  stream 
surface  is  totally  cylindrical  shape,  the  governing  equation 
to  be  used  for  the  cascade  analysis  will  be  a  Laplace 
equation.  Unfortunately,  the  stream  surface  is  of  three- 
dimensional  cone  shape  in  general  for  the  tail  cone  section 
of  the  underwater  vehicle.  The  field  governing  equation  now 
becomes  a  Poisson  equation,  and  the  results  of  powerful 
potential  theory  analysis  are  no  more  applicable.  A  method 
of  correcting  the  effect  of  the  Poisson  equation  on  the 
potential  theory  results  is  introduced  to  modify  the  blade 
profile  shape  obtained  in  the  potential  theory.  In  choosing 
the  blade  profile  shape,  the  experimental  data  are  used  to 
ensure  that  there  is  no  chance  of  flow  separation  due  to 
overloading  on  the  blade.  Furthermore,  based  on  the  calcu¬ 
lated  velocity  along  the  blade,  the  cavitation  inception  is 
checked.  If  there  exists  a  chance  of  either  flow  separation 
or  cavitation,  the  loading  distribution  from  hub  to  tip 
should  be  changed.  If  such  a  change  is  made,  and/or 
thickness  of  blades  is  changed,  the  streamline  curvature 
method  should  be  used  again  to  determine  the  new  location  of 
streamline  or  stream  surface.  This  iterative  procedure  is 
to  be  repeated  until  an  overall  convergent  solution  is 
obtained.  Section  3.3  describes  the  technical  approach  to 
be  used  for  the  blade-to-blade  flow  analysis.  The  two- 
dimensional  theory  applicable  to  the  mapped  stream  surface 
is  summarized  in  Section  3.2.  The  overall  product  of  three- 
dimensional  pumpjet  design  is  described  in  Section  3.4. 


3.1 


BLAOE-THROUGH  FLOW  ANALYSIS  -  STREAMLINE  CURVATURE 
METHOD 


The  mathematical  formulation  for  the  streamline  curvature 
method  (SCM)  has  been  described  in  the  1985  report  (Furuya 
and  Chiang,  1986). 

3.2  TWO-DIMENSIONAL  ANALYSIS 

In  the  following,  the  major  concept  of  Mellor’s  two- 
dimensional  cascade  theory  will  be  summarized  for  the 
reader’s  convenience. 

3.2.1  Linearized  Cascade  Theory 

A  two-dimensional  cascade  theory  is  implemented  here  to  be 
used  for  calculating  the  lift  coefficient  for  the  three- 
dimensional  flow  discussed  in  Section  3.3.  The  boundary 
conditon  (Eqn.  3.2-14)  has  to  be  modified,  as  discussed  in 
Section  3.3,  when  this  cascade  theory  is  applied  in  the 
quasi  three-dimensional  flow. 

The  lift  coefficient  is  determined  for  any  given  cascade 
geometry  which  is  specified  by  the  solidity  (blade  chord-gap 
ratio,  c/s)  and  the  stagger  angle,  X  (Figure  3-2).  Symbols 
used  in  Mellor  (1959)  were  followed;  use  Li  and  C^  in 
denoting  the  ideal  lift  force  and  lift  coefficient  when  the 
drag  is  zero.  Then  (see  also  Weinig,  1964) 

Li  =  ( s/c ) pWmAV^  (3.2—1) 

and  CL.  =  (Li/c)/(rPwm2) 

=  2(s/c)(AVe/Wm)  (3.2-2) 

where  s  is  the  blade  pitch,  c  is  the  chord  length,  p  is  the 
fluid  density,  Wm  is  the  mean  relative  velocity,  and  AVe  is 
the  difference  between  the  peripheral  velocity  at  the  exit 
and  that  at  the  inlet. 


Replacing  the  product  sAVe  by  the  line  integral  fv-dr  on  a 
closed  path  comprising  two  streamlines  s  distance  apart  and 
joined  by  two  lines  parallel  to  the  •e-direction,  we  have 
(Weinig,  1964)  . 


where 


cL.  =  2r/(cwm) 


I*  ■  f>v  •  dr 


(3.2-3) 

(3.2-4) 


is  the  circulation  around  a  profile  (Wislicenus,  1965). 

The  camber  is  assumed  to  be  sufficiently  small  so  that  the 
chord  length  is  substantially  equal  to  the  distance  measured 
along  the  camber  line.  Then,  the  circulation  around  a  thin 
wing  profile  is  given  by  (Abbott  and  von  Doenhoff ,  1959) 


T  =  Sl  y  dx 


(3.2-5) 


where  y  is  the  difference  in  velocity  between  the  suction 
and  pressure  surfaces,  which  is  also  the  strength  of  the 
vortex  sheet  comprising  the  blade  camber  line  (von  K^rm^n 
and  Burgers,  1963).  Therefore,  equation  (3.2-3)  is 
expressed  by 


cl_i  =  2  /J  (y/Wm)  d  (x/c) 


(3.2-6) 


The  cambered  blade  is  built  up  by  superimposing  vortices  on 
the  camber  line  and  a  distribution  of  sources  and  sinks  on 
the  camber  line  to  account  for  the  profile  thickness 
effects.  The  distribution  of  source  (sink),  q,  is  (Mellor, 
1959) 


q  =  wm  dyt/dx  +  d(uyt)/dx 


(3.2-7) 


where  the  thickness  of  blade  is  denoted  by  yt  and  the 
induced  chord-wise  velocity,  u,  is  considered  constant  along 
the  y-direction  within  the  profile.  The  second  term  can  be 
shown  to  be  negligible  (Mellor,  1959)  and  so  we  have 


=  (t/c)ft(x/c)  (3.2-8) 

J 

where  ffc(x)  ■  3ffc/3x  (3.2-9) 

and  the  thickness  function  fb  is  defined  by 

yt/c  ■  (t/c)ft(x/c)  (3.2-10) 

where  t  is  the  maximum  thickness  of  the  blade. 

The  camber  function  fc  is  defined  by 

yc/c  ■  Cb  fc(x/c)  (3.2-11) 

where  yc  denotes  the  camber  distribution  and  Cb  is  defined 
by 


Cb  ■  2  /J  (dyc/dx)cose  de  (3.2-12) 

in  which 

cose  ■  1  -  2  x/c  (3.2-13) 

A  blade  is  approached  by  a  mean  velocity  Wm  at  a  mean  angle 
am.  To  satisfy  the  condition  that  the  normal  velocity 
vanishes  at  the  boundary,  the  flow  velocity,  together  with 
the  induced  velocity,  should  be  tangent  to  the  surfaces. 
Neglecting  the  thickness  effect,  the  boundary  condition  at 
x0  becomes 

(Wm  sin  am  +  v0 )/( Wm  cos  a,*  +  u0)  =  (dyc/dx)0 

*  Cb  f'Q  (xQ/c)  (3.2-14) 

where  u0  and  v0  denote  respectively  the  x-  and  y-components 
of  the  induced  velocity  at  x0  on  the  O^h  blade  with  x 
measured  along  the  chord  from  the  leading  edge. 


To  find  the  lift  coefficient  by  equation  (3.2-6),  we  assume 
that  y/Wm  may  be  represented  by  a  trignometric  series  as 
(Abbott  and  von  Doenhoff,  1959) 


y/Wm  =  2A  (1+cos  e)/sin  «  +  4  2-An  sin  n-e  (3.2-15) 

which  is  zero  at  the  trailing  edge  of  e  =  tt  so  that  the 
Kutta  condition  is  satisfied.  This  distribution  of  vor¬ 
tices,  together  with  the  distribution  of  sources/sinks  shown 
in  equation  (3.2-8),  may  be  used  to  obtain  the  components  of 
induced  velocity,  u0  and  v0.  The  components  of  induced 
velocity  are  then  substituted  into  equation  (3.2-14).  With 
the  aid  of  (von  Kerman  and  Burgers,  1963;  Milne-Thomson , 
1966) 

S0  [cos  n«/(cos  «■  -  cos  d0)J  d-©-  =  ir  sin  n-»0/sin 

(n  =  0,  1,  2,  . . . )  (3.2-16) 

the  following  equation  is  obtained: 

00  , 

Eq  Angn  =  sinam  -  cb  fc(s„)  cosam 

*  Cb  £  Anhn  -  |  (CbB-T)  (3.2-17) 

n=o 

where  gn,  hn,  B,  and  T  are  defined  in  Mellor  (1959)  except 
that  ft(0o)  should  be  replaced  by  f{(s)  in  defining  T.  An 
are  the  Fourier  coefficients  to  be  evaluated. 

For  a  cascade  of  certain  solidity,  incident  angle,  and 
stagger  angle,  equation  (3.2-17)  is  used  to  calculate  N 
coefficients,  An,  based  on  the  camber  and  thickness  at  N 
locations  along  the  chord.  This  study  follows  the  solution 
method  of  Mellor  (1959)  which  greatly  reduces  the  calcula¬ 
tion  labor  when  a  set  of  solutions  as  functions  of  the 
solidity  and  stagger  angle  is  desired. 

Equation  (3.2-17)  is  multiplied  by  cos  ke0  and  integrated 
from  0  to  ir  to  obtain  an  equation  for  numerical  integration. 

Having  the  double  Fourier  integral  functions  calculated,  the 
cascade  coefficients,  An,  were  computed.  Then  the  lift 
coefficient  was  obtained  from 


C|_i  =  2it(Aq  +  Ai) 


(3. 2-18) 


which  is  derived  by  inserting  equation  (3.2-15)  to  equation 
(3.2-6) . 


3.2.2  Data  Analysis 

The  NACA  65-series  experimental  data  given  by  Herrig,  et  al. 
(1951)  are  used  to  do  the  multiple  regression  analysis  in 
the  present  study.  The  data  include  the  cascade  lift  coef¬ 
ficient  and  the  design  angle  of  attack.  The  design  angle  of 
attack  is  a  function  of  solidity  and  blade  camber  (Figure 
3-3).  The  lift  coefficient  is  a  function  of  stagger  angle, 
solidity,  blade  camber,  and  angle  of  attack.  Among  the  lift 
coefficient  values,  only  those  associated  with  the  design 
angle  of  attack  are  used  in  the  data  analysis. 

Based  on  28  data  points  in  Figure  3-3,  the  fourth  order 
polynomial  equation  of  the  design  angle  of  attack,  ad, 
obtained  by  the  multiple  regression  analysis  is 

ad  =  -  1.78681  a4  +  0.51975  a3Cb  +  0.02078  a2Cb2 

-  0.06408  aCb3  +  0.01716  Cb4  +  7.68063  a3 

-  2.50213  a2Cb  +  0.30280  aCb2  +  0.01298  Cb3 

-  12.95580  a2  +  6.54778  aCb  -  0.30496  Cb2 
+  12.87888  a  +  2.78086  Cb  -  2.22656 

(3.2-19) 

where  a  denotes  solidity  and  Cb  represents  the  camber. 
Figure  3-4  shows  the  correlation  between  the  calculated 
value  from  Eqn.  (3.2-19)  and  the  experimental  data.  The 
mean  residual  is  -0.00002,  the  standard  deviation  is  0.04, 
and  the  maximum  residual  is  0.09  degrees  which  is  smaller 
than  the  error  bound  of  the  original  data. 

Cascade  lift  coefficient  at  the  design  angle  of  attack  has 
79  data  as  a  function  of  stagger  angle,  solidity,  and  camber 
in  Herrig,  et  al.  (1951).  Two  data  at  the  falling  limb,  in 
the  figure  of  lift  coefficient  vs.  angle  of  attack,  are 


removed  from  the  sample.  A  total  of  77  data  is  used  in  the 
multiple  regression  analysis.  The  resultant  fourth  order 
polynomial  equation  is 

Cb  =  0.00000038594  B4  -  1.66604  a4  -  2.74752  C £4 

-  0.000050739  8 3a  -  4.42872  a3 Ct  -  0.093633  C^3B 

+  0.0013187  02a2  -  10.87314  a2C^2  -  0.0016033  Ct2 B2 

-  0.0014369  8a3  -  9.83274  aC^3  -  0.000035308  C^B3 

-  0.000020018  83  +  8.26201  a3  +  20.91537  C^3 

+  0.0040173  82c  +  26.19466  a2C^  +  0.33569  C £28 

-  0.11946  Ba2  +  41.23668  oC^2  +  0.0059786  C^B2 

-  0.00080978  B2  -  14.15108  a2  -  45.84094  ZL2 

-  0.010265  Bcr  -  47.03403  oC^  -  0.39504  C^B 
+  0.015917  B  +  13.59117  a  +  36.29869  C* 


-  5.23047 


(3.2-20) 


The  results  calculated  from  this  equation  are  shown  as  solid 
lines  in  Figures  3-5  to  3-8  for  different  relative  flow 
angle  at  the  blade  inlet.  Also  shown  in  these  figures  are 
discrete  data  at  different  conditions.  The  results  from  the 
regression  analysis  fit  well  with  the  data. 

3.3  THREE-DIMENSIONAL  ANALYSIS  - 
BLADE-TO-BLADE  FLOW 

This  section  illustrates  the  theory  and  procedure  to  solve 
the  blade-to-blade  flow  on  each  stream  surface  in  an  axisym- 
metric  three-dimensional  flow  environment. 


3.3.1 


Differential  Equations 


Under  the  assumption  that  an  axisymmetric  stream  surface 
exists  in  a  rotating  machine,  from  the  conservation  equation 
of  steady  circulation,  i.e,  7  x  w  +  2o>  =  0,  the  following 
relation  is  obtained  for  the  relative  flow, 


3(rwQ)  3r 

—t— - =  2o) r  r—  =  2ojr  sinX 


(3.3-1) 


where  wm  and  we  are  relative  flow  velocities  in  the  direc- 

» 

tions  of  m  and  e,  r  measures  the  radial  distance,  and  \  is 
the  angle  of  the  line  tangent  to  the  stream  surface  at  the 
point  of  interest  made  with  the  axis  of  rotation  (Figure 
3-9).  The  continuity  equation  for  the  same  stream  surface 
is  also  written 


3(bpwQ)  3(bprwm) 
3e  +  3m 


(3.3-2) 


where  b  is  the  thickness  of  stream  surface  and  p  is  the 
fluid  density. 


Then,  a  stream  function  can  be  defined  by 


W  ;  ri  21 

9  Bp  3m  ’ 


_i  ajJ_ 

bp  r3e 


(3.3-3) 


Substitution  of  w9  and  wm  in  Eqn.  (3.3-3)  into  Eqn. 
(3.3-1)  yields 


+  F  ’  B  Im  fm  =  ”2bpu)  sinX 


(3.3-4) 


3.3.2  Transformation 


Consider  a  Cartesian  coordinate  system  (X,Y)  with  the  origin 
0  at  the  leading  edge  of  a  blade  and  the  X-axis  in  the  axial 
direction  (Figure  3-2),  the  three-dimensional  axisymmetric 
stream  surface  given  by  Eqn.  (3.3-4)  can  be  mapped  onto  this 
two-dimensional  X-Y  plane  by  the  following  mapping  functions 

dx  ro  dY  f  0 
=  —  ,  rnr  =  -r„  (3.3-5) 


where  r0  is  an  arbitrary  constant  which  is  used  for  the  pur¬ 
pose  of  scaling  between  the  physical  coordinate  space  and 
mapped  plane  (X,Y).  By  using  Eqn.  (3.3-5),  the  governing 
equation  (3.4-4)  can  now  be  written  in  the  (X,Y)  coordinate 
system  as 


=  -2bp0J  C  — )  sinX 
o 


.  l  Mbp)  a£  3(bp)  s£ 
Bp)  3X  3X  3  Y  3  Y  l  . 


(3.3-6) 


Also,  the  relative  velocities  in  the  X-  and  Y-directions  are 
given  by 


fa)  w  -  *  It  -  £_  w  -  £-  c 

Lajwx  Ep3Y  r  m  ’  r  Sn 

oo 

(b)  wY  =  -  -  f-  w0  =  J_  (c0  -u) 


(3.3-7) 


where  cm  and  ce  are  absolute  flow  velocities  along  m  and  e 
directions,  respectively. 

On  the  mapped  plane,  the  relative  flow  angles  at  the  inlet 
and  exit  of  a  blade,  8i  and  82*  respectively,  are  obtained 
from  the  equivalent  velocity  diagram  to  be 


tan  ei  -  F 


(3.3-8) 


tan  82  *  F 


(3.3-9) 


where  the  subscripts  1  and  2  denote  the  condition  at  the 
inlet  and  exit,  respectively,  of  a  blade,  and  wx®  is  the 
mean  value  of  wxi  and  wx2- 

3.3.3  Effects  of  Streamline  Inclination  and 
Meridian  Velocity  Variation 

As  seen  from  Eqn.  3.3-6,  the  governing  equation  for  the 
( X , Y )  plane  is  now  a  Poisson  equation  instead  of  the  Laplace 
equation  which  exists  only  for  a  flow  on  a  perfectly  cylin¬ 
drical  stream  surface  with  uniform  velocity  distribution. 
Therefore,  the  results  obtained  from  the  two-dimensional 
linear  cacade  theory  should  be  corrected  according  to  the 
right-hand  side  term  of  Eqn.  (3.3-6).  It  is  readily 


understood  that  these  right-hand  side  terms  are  satisfied  by 
distributing  the  following  vortices,  £ ,  and  sources,  y,  on 
the  entire  (X,Y)  plane 


r  ^ 

(a)  c,  =  (7  x  w)x  y  =  2u)  (— )  sinX 


(3 . 3t10 ) 


( b)  y  =  (7 • w) x  y 


,  .1.  htbp)  it  3(bp)  a£ 

[bp}-1  )  3X  3 Y  3 Y  3X 


By  adding  the  induced  velocities  calculated  from  t,  and  y, 

the  blade  profile  shape  or  equivalently  the  camber  obtained 

in  the  conventional  two-dimensional  analysis  will  be 

corrected.  It  should  be  noted  that  the  first  term  on  the 

} 

right-hand  side  of  Eqn.  (3.3-6)  arises  from  non-zero  X  , 
i.e.,  the  stream  surface  is  not  parallel  to  the  axis  of 
rotation,  whereas  the  second  group  of  terms  is  due  to  the 
non-uniform  thickness  of  stream  surface  or  tube  caused  by 
the  variation  of  meridional  velocity.  Needless  to  say,  if 

9 

X  =  0  and  bp  is  constant,  Eqn.  (3.3-6)  becomes  a  Laplace 
equation  and  thus  a  two-dimensional  linear  cascade  theory 
holds . 

A  method  similar  to  the  present  one  was  developed  by  Inoue 
and  his  colleague.  In  their  study  (e.g.,  Inoue,  et  al., 
1980),  there  exist  a  few  major  drawbacks,  some  of  which 
could  potentially  lead  to  a  substantial  error  in  the  final 
design.  First  of  all,  since  they  use  a  two-dimensional 
linearized  cascade  theory,  the  error  becomes  significant  for 
high  solidity  and  high  stagger  angle  area,  i.e.,  near  the 
hub,  although  they  introduce  experimental  data  in  a  later 
step  of  the  analysis.  Secondly,  their  velocity  triangle 
used  for  determining  the  incoming  flow  angle  to  the  blade  is 
in  error  of  the  first  order  since  they  did  not  take  into 
consideration  the  effect  of  non-cy lindrical  and  variable 
thickness  stream  surface.  Finally,  due  to  the  use  of  the 
linearized  cascade  theory,  they  failed  to  obtain  the  veloc¬ 
ity  distribution  and  therefore  a  boundary  layer  analysis 


i.ui  up  m  f  nwwmiy  mi'iin  uwin'jn  w^i^!^wnuiAMi''rvft  *  *»  w  ■  ■  ■■  ■  «■■■■■■  -m  ■  v 
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and  cavitation  inception  analysis  are  not  possible. 

With  these  aspects  in  mind,  effort  has  been  made  in  the 
current  GHR  project  to  improve  the  accuracy  of  the  linear 
cascade  theory  as  well  as  to  avoid  the  singular  behavior  of 
velocity  at  the  leading  edge  of  blade.  Detailed  discussions 
on  the  loading  correction  and  leading  edge  correction  have 
been  presented  in  the  FY  ’85  Report. 

3.3.4  Induced  Velocities 

If  the  inclination  of  stream  surface  is  small  such  as  that 
in  an  axial-flow  case,  an  approximation  solution  of  the 
velocity  induced  by  the  distributed  vortices  (Eqn.  3.3-10a) 


are  obtained  by 
et  al.,  1979): 

a 

replaced  average 

vorticity  £  (Inoue, 

1 

f 

COSX 

£dx 

C 

c  cosX 

J 

0 

uo 

r  2- 
r2 

(3.3-11) 

c  cosX 

r 

2 

0 

where  c 

is  the  chord 

length,  u0  is  the 

speed  of  blade  at  the 

reference  radius  r0,  and  subscripts  1  and  2  indicate  the 
inlet  and  exit,  respectively,  of  the  blade. 

Similarly,  when  the  variation  of  axial  velocity  is  small, 
the  distribution  of  sources  (Eqn.  3.3-10b)  is  replaced  by  a 
uniform  distribution: 


-  _  WX2~WX1 
u  “  c  cosX 

-  A"x 
c  cosX 

where  w  denotes  the  average  relative  velocity. 


(3.3-12) 
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.y  '  /.y>: 


v-  j-. 


clvl'rfiVl'jCvCv>lvl'w>lvlv 


'•(r'vjl-lvlv: 


Consider  another  Cartesian  coordinate  system  (x,y)  with  the 
origin  at  0  and  the  x-axis  in  the  chordwise  direction,  which 
has  a  stagger  angle  X  relative  to  the  axial  direction 
(Figure  3-4).  The  mean  flow  velocity  along  the  chord  direc¬ 
tion  is 


cosX  +  wVm  sinX 
X®  X®  Y® 


=  wXoo  cosX  (1  +  tang^  tanX).  (3.3-13) 

The  induced  velocities,  relative  to  wx®,  due  to  the  uniform 
distribution  of  vortices,  ?,  and  sources,  u,  are 


where 


x® 


V 

w 

X® 


y  r*  _  1  tan  X)  _ tan  X _ 

x  Lc  c  J  1  +  tanX  tang 


x  tf  -  I  tan  x>  1  +  tan 


tanX  tang 


IM  _  r  r*  .  I  tan  X)  _ L 

wxoo  4  L  c  c  1  +  tanX 


tang. 


5“  *  -  \  »  — 

X® 


tan  X 


tanX  tang 


C  c  cosX 
"x® 


(3.3-14) 


(3.3-15) 


(3.3-16) 


(3.3-17) 


_  2  _  2 
1  2  “1 


(3.3-18) 
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m  i 


r„  Uo 
0  0 


(3.3-19) 


are  streamline  inclination  parameter  and  axial  velocity 
variation  parameter,  respectively,  and 
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rl  "ml 


"  r2  "m2 


2  ro 

o  o 


(3.3-20) 


is  a  local  flow  coefficient,  with  subscript  m  denoting  the 
meridional  component  of  the  velocity. 


Eqns .  (3.3-14)  and  (3.3-15)  are  good  only  if  the  streamline 
inclination  is  small.  Eqns.  (3.3-16)  and  (3.3-17)  are 
obtained  by  ignoring  the  blockage  effect  of  blade  thickness. 
In  the  following  discussion,  the  blockage  effect  is  con¬ 
sidered  and  the  induced  velocity  from  the  distributed  vor¬ 
tices  given  by  Eqn.  (3.3-  10a )  are  obtained  by  solving  the 
Poisson  equation 


V2lp  =  -2bpoj  (  y~)  sin* 
o 


to  give  the  solution 


VCY  *  U0 


0 


2  2 
2  rl  *  r 2 

'  2r  2 


where  u0  is  a  reference  velocity. 


(3.3-21) 


(3.3-22) 


The  induced  velocity  due  to  the  distribution  of  sources 
given  by  Eqn.  (3.3-lQb)  is  approximated  by  (Inoue,  et  al  .  , 
1980) 


j-  .  t 


-T- 
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(3. 3-23 ) 


with  the  blockage  factor  of  blade  thickness,  ,  put  into 
the  consideration. 


By  decomposing  both  vrY  and  v^x  into  the  x  and  y  directions, 
Eqns.  (3.3-22)  and  (3.3-23),  together  with  (3.3-13),  become 


/  tanX 


(3.3-25) 


Flow  Skewness  in  Diagonal  Contracting  Channel 


The  major  goal  of  the  proposed  theory  is  to  make  corrections 
of  flow  stream  tube  nonuniformness  and  diagonal  flow  effects 
on  the  two-dimensional  axial  cascade  flow.  This  is  per¬ 
formed  within  the  framework  of  perturbation  method  by  uni¬ 
formly  distributing  vortices,  y,  and  source/sink ,  c,  ,  '  in 
cascade  row,  as  described  in  Section  3.3.4.  The  blade  pro¬ 
file  shape  will  be  properly  modified  in  order  to  take  y  and 
T,  into  account.  Due  to  the  distribution  of  uniform  vortices 
and  source/sink,  the  flow  skewness  problem  exists,  i.e.,  the 
first  order  effect  on  the  flow  incidence  angle  is  generated 
by  such  singularity  distribution. 

As  shown  in  Figure  3-10,  the  constant  vortices,  £,  are 
distributed  in  a  strip  over  the  entire  cascade  row.  The 
complex  potential  W  for  a  point  vortex  c  placed  at  z'  is 
expressed 


W  =  tn  (z-z’  ) 


(3.3-293 


and  the  induced  velocities  are  obtained  by  taking  the 
deviations  of  W 


u^  -  iv5 


dW  i_C  x-x*  +  iy* 
^  =  2*  (x-x’  F+y’’ 


(3.3-30) 


The  induced  velocities  due  to  the  uniform  distribution  of 
vortices  over  a  strip  of  cascade  region  are  obtained  by 
integration 


lcos\  ® 

\  ~  \  J  dx’  dy1 

0  -00 


tcosX  00  —  „  > 
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—  £cosX 

T"  /  (-  ?-0)  dx’ 


=  ^  £cosX , 


(3.3-31) 


whereas  x+  +® 


v  =  -  -j  £cosX.  (3.3-32) 

Therefore,  by  distributing  If  uniformly  as  described  above, 
the  upstream  and  downstream  flows  are  equally  deviated  from 
the  geometric  mean  flow,  w«  (see  Figure  3-10  again). 

On  the  other  hand,  a  more  rigorous  mathematical  model  by 
Mani  and  Acosta  (1968)  will  lead  to  a  totally  different 
conclusion.  Since  this  point  has  been  ignored  by  many 
researchers,  let  us  review  herein  the  details  of  Mani  and 
Acosta’s  work.  In  their  work,  the  channel  contraction  was 
chosen  to  be  of  an  exponential  shape  as  shown  in  Figure 
3-11.  By  satisfying  the  proper  boundary  conditions,  the 
velocity  component  in  the  y-direction  due  to  a  point  vortex 
of  unit  strength  at  the  origin  is  obtained  (see  the  report 
of  Mani  (1966)), 
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for  x’  <-b 
2,.  >  2 


(nfr.)  »  «  tn  (3.3-3*) 

+  0  (a^)  for  x’  >a  . 


Therefore,  the  induced  velocity  due  to  a  distribution  of  the 
vortices  over  a  chord  length  in  the  cascade  configuration  is 
obtained 
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(3.3-39) 


The  second  terms  of  (3.3-36)  and  (3.3-37)  are  defined, 
respectively 
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direction.  If  E  =  2a  =  2b, 
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(3.3-40) 


Rewriting  this  induced  velocity  in  two  parts,  by  vp  and  va, 


v  = 


vr  "  Va ’  as  x'*‘  "® 


-v  -  v  ,  as  x+  +oo 
P  o 


where 


The  velocity  diagram  is  now  written  in  Figure  3-12.  It  is 
clearly  seen  that  there  exists  an  error  of  first  order  of  a 
if  woo  is  chosen  as  the  reference  velocity.  Instead,  w<^ 
should  be  chosen  as  the  reference  velocity  so  that  the  blade 
setting  will  be  changed.  This  error  is  substantial  since 
the  three-dimensional,  diagonal  flow  correction  handled 
under  the  current  method  is  order  of  magnitude  a.  This  fact 
indicates  that  not  only  the  blade  camber  profile  needs 
correction  as  discussed  earlier  but  also  the  blade  setting 
should  be  changed  in  order  to  take  such  3-D  correction  into 
account  correctly.  Inoue,  et  al.  (1979,  1980)  did  not  see 
this  point  in  their  series  of  papers. 


There  exist  two  possible  ways  of  implementing  this  idea: 

1)  by  assessing  the  flow  channel  contraction  or  expansion  in 
terms  of  equivalent  a,  Equation  3.3-40  will  be  used,  and 

2)  a  more  accurate  calculation  will  be  made  by  using  the 
distributed  vortices  over  the  cascade  strip.  The  method  of 
actually  incorporating  this  flow  "skewness”  into  the  current 
design  theory  will  be  one  of  the  major  tasks  in  the  FY-87 
GHR  project. 


3.3.7  Secondary  Flow  Correction 

The  secondary  flow  theory  commonly  used  for  axial  flow  is 
applicable  to  the  current  diagonal  flow.  The  strength  of 
vortex  at  the  exit  of  blade  channel  due  to  the  secondary 
flow,  Csf»  expressed 

N  [r  ^1  I  ds  _1 dT 

^  SF  2irr  2cos8  2  r1  W1  b  ^  ”  #  ”  cose  dq~ 

where  N  =  number  of  blades, 

; .  =  the  component  of  vortex  normal  to  the  direction 
of  relative  flow  on  the  rotating  surface, 

w^  =  relative  incoming  velocity 

s  =  line  element  along  the  blade  on  the  rotating 
surface 

w  =  relative  flow  velocity  on  the  blade  surface 
r  =  circulation  of  blade 

q  =  q-line  used  for  the  streamline  curvature  method 

e  =  angle  of  a  tangent  line  of  meridional  line  made 
with  q-line 

b i ,  b2  =  widths  of  flow  at  inlet  and  exit,  respectively 

82  =  flow  angle  at  exit  relative  to  chord  line. 

By  distributing  the  above  secondary  vortex  at  the  exit  of 
blade  channel  and  satisfying  the  boundary  condition  on  the 
hub  and  casing,  the  induced  flow  effect  on  Ce  (called  ACe 
hereafter),  will  be  calculated  by  solving  the  corresponding 
Poisson’s  equation.  In  the  current  design  problem,  this  ACe 
will  be  taken  into  consideration  for  determining  the  blade 
profile  shape.  More  detailed  numerical  analysis  will  be 
developed  in  future  work. 

3.4  DESIGN  PROCEDURE 


The  design  procedure  presented  here  is  based  on  the  assump¬ 
tion  that  the  shape  and  steady  velocity  of  underwater 


vehicle  is  given.  The  whole  procedure  can  be  repeated  for  a 
modified  shape  or  velocity.  At  the  beginning,  some  parame¬ 
ters,  such  as  the  location  of  the  rotor,  the  number  of 
blades,  and  rotation  speed,  have  to  be  chosen  based  on 
experiences . 

Then  [Figure  3-1),  the  stream  curvature  method  (SCM)  is  used 
to  determine  the  streamlines  for  a  through-flow  in  a  merid¬ 
ional  plane.  As  discussed  in  Section  3.1,  the  method  is 
based  on  solving  the  momentum  equation  along  quasi- 
orthogonal  lines  (q-lines)  and  the  conservation  of  mass  is 
kept  along  the  meridional  direction. 

After  streamlines  are  determined,  average  stream  surfaces 
are  taken  to  be  the  revolution  of  streamlines  about  the  axis 
of  rotation.  Then  a  program,  DSN3,  is  used  to  do  the 
remaining  design  procedure.  Figure  3-13  depicts  the  macro 
view  of  the  procedure  while  Figures  3-14  through  3-18  show 
details  of  each  substep  based  on  the  theory  presented  in 
Sections  3.1  to  3.3. 

Figure  3-14  shows  the  input  data  required  in  the  general 
background.  Input  data  related  to  each  individual  cross- 
section  are  to  be  read  in  Subroutine  INP2  as  shown  in  Figure 
3-15.  Some  input  data  are  obtained  from  the  meridional  flow 
solutions  computed  in  the  SCM  program. 

Figure  3-16  shows  the  subroutine  to  calculate  the  flow  vel¬ 
ocities  in  the  mapped  (X,Y)  plane  as  discussed  in  Section 
3.3.2.  The  subscripts  1  and  2  denote  conditions  at  inlet 
and  exit,  respectively,  of  a  blade. 


One  basic  concept  in  the  present  design  procedure  is  that 
the  blade  camber  design  based  on  cascade  experimental  data 
is  adjusted,  by  considering  the  effects  of  streamline  incli¬ 
nation  ana  nonuniform  velocity  distribution,  such  that  both 
the  final  flow  turning  angle  and  the  incident  angle  relative 
to  the  axis  of  rotation  are  the  same  as  those  obtained  in 
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the  two-dimensional  flow  without  the  three-dimensional  fac¬ 
tors. 

The  camber  and  turning  angle  for  the  case  of  uniform, 
parallel  flow  are  evaluated  based  on  the  required  lift  coef¬ 
ficient  and  solidity  (Figure  3-17).  The  method  relies  on 
the  experimental  data  presented  in  Section  3.2.2. 

In  the  presentation  of  three-dimensional  effects,  which  are 
to  be  considered  in  the  boundary  condition  as  discussed  in 
Section  3.3.5,  the  camber  and  stagger  angle  are  adjusted 
until  the  turning  angle  in  the  mapped  two-dimensional  plane 
is  the  same  as  that  with  original  camber  in  the  uniform, 
parallel  flow  (Figure  3-18).  A  linearized  cascade  theory 
(Section  3.2.1),  together  with  proper  boundary  condition 
which  has  the  three-dimensional  effects  considered  (Section 
3.3.5),  is  used  to  calculate  the  cascade  lift  coefficient 
and  the  associated  total  circulation. 

If  the  calculated  turning  angle  for  the  three-dimensional 
flow  case  is  different  from  the  desired  value,  the  camber 
and  stagger  angles  are  adjusted  in  an  iteration  process 
until  the  desired  value  is  achieved  to  within  certain 
tolerance  criterion.  When  the  result  is  converged,  the 
final  camber  is  the  one  which,  under  the  influence  of  three- 
dimensional  flow,  will  yield  the  desired  lift  coefficient. 


4.0 


CONCLUSIONS 


During  the  FY-86  GHR  study,  the  mathematical  model  of  the 
blade-to-blade  theory  for  the  three-dimensional  pumpjet  flow 
has  been  developed.  The  stream  surface  of  nonuniform 
thickness  in  the  diagonal  flow  was  mapped  into  plane  by  a 
functional  mapping  function.  The  rotor  blades  now  become  a 
row  of  blades,  i.e.,  cascade  configuration.  The  Laplace 
equation  in  the  three-dimensional  coordinate  system  is  then 
converted  into  a  Poisson’s  equation.  This  Poisson’s 
equation  is  solved  as  the  corresponding  Laplace  equation 
combined  with  the  source/sink  and  vortex  distribution  on  the 
rotating  flow  field  as  a  correction.  More  specifically,  the 
two-dimensional  cascade  solution  is  corrected  in  terms  of 
blade  camber  profile  shape  due  to  the  induced  velocities 
generated  by  the  source/sink  and  vortex  distribution.  An 
iteration  procedure  is  necessary  to  satisfy  the  boundary 
condition  on  each  blade,  namely,  the  flow  tangency  con¬ 
dition  since  the  original  two-dimensional  cascade  solution 
will  not  satisfy  the  tangency  condition  after  the  correction 
is  made. 

An  additional  important  correction  due  to  the  correction  by 
the  singularity  distribution,  particularly  by  the  vortex 
distribution,  is  that  of  flow  skewness.  It  seems  that  this 
point  was  ignored  by  most  researchers,  but  is  as  important 
as  the  camber  correction.  It  is  of  first  order  of  three 
dimensionality.  This  skewness  changes  the  flow  directions 
at  upstream  and  downstream  in  the  same  direction  so  that  the 
blade  setting  should  also  be  changed.  It  means  that  the 
stagger  angle  of  the  cascade  blade  should  be  changed  prop¬ 
erly. 

During  the  FY-87  GHR  study,  these  two  theories,  i.e.,  blade- 
through  flow  theory  (i.e.,  Streamline  Curvature  Method)  and 
blade-to-blade  flow  theory,  will  be  combined  and  the  itera¬ 
tion  procedure  will  be  established.  After  developing  com¬ 
puter  codes,  a  sample  pumpjet  design  will  be  attempted. 
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FIGURE  1-3  A  typical  load  distribution 
in  terms  of  Vc  for  pumpjet 
rotor  blade  where  V,  is  cir 
cumferential  component  of 
the  turned  flow  velocity 


FIGURE  1-4  Typical  pump jet  rotor  blade 
figuration,  (a)  top  view  ani 
(b)  upstream  view 
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FIGURE  3-1  Floy*  Chart  of  the  Selected 
Pumpjet  Design  Method 
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4  Comparison  of  design  angle  of  attack 
obtained  from  the  multiple  regression 
analysis  and  those  obtained  from  the 
laboratory 
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FIGURE  3- 


5  Camber  as  a  function  of  cascade  lift 
coefficient  and  solidity  obtained  from 
regression  analysis  data  (solid  line) 
compared  with  original  data  (discrete 
data  point) ,  for  31  =  30° 
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FIGURE  3-6  Camber  as  a  function  of  cascade  lift 

coefficient  and  solidity  obtained  from 
regression  analysis  data  (solid  line) 
compared  with  original  data  (discrete 
data  point),  for  3^  =  45° 
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FIGURE  3- 


7  Camber  as  a  function  of  cascade  lift 
coefficient  and  solidity  obtained  from 
regression  analysis  data  (solid  line) 
compared  with  original  data  (discrete 
data  point) ,  for  6,  =  60° 
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Figure  3-13 


General  flow  chart  to  design 
blade  in  a  flow  of  three- 
dimensional  character 
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Figure  3-14  Flow  chart  for  Subroutine  INP 


SUBROUTINE  INP2' 
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Figure  3-15  Flow  chart  for  Subroutine  INP2 


Figure  3-16  Flow  chart  for  Subroutine  MAP32 
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Figure  3-17  Flow  chart  for  Subroutine  DSN2 
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Figure  3-18  Flow  chart  for  Subroutine  PFM3 


